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ABSTRACT 


The  low  Reynolds  number  hypersonic  flight  regime  introduces  the 
aerothermochemical  interaction  of  low  density  hypersonic  flow  and 
real  surfaces.  For  the  ballistic  re-entry  of  slender  bodies  or  for 
sustained  hypersonic  flight  at  high  altitudes,  it  may  prove  advantagous 
to  use  local  mass  transfer  cooling  over  certain  portions  of  the  sur¬ 
face  to  reduce  the  surface  temperature  where  it  is  expected  that 
oxidation  of  the  surface  will  otherwise  produce  excessive  shape  changes. 
The  effusion  of  air  through  a  porous  surface  has  been  considered  in 
this  investigation,  and  the  simultaneous  heat  and  mass  transfer  proc¬ 
esses  experienced  by  hypersonic  vehicles  flying  at  high  altitudes  have 
been  analyzed.  Because  stagnation  region  flows  represent  the  sever¬ 
est  heating  problems,  numerical  results  are  presented  which  are  ap¬ 
plicable  at  the  forward  stagnation  point  of  axi- symmetric  hypersonic 
vehicles. 

The  governing  equations  utilized  in  the  present  analysis  are  the 
equations  of  change  for  the  flow  of  a  compressible  chemically  reacting 
gas.  Included  were  the  conservation  equations  of  mass,  momentum 
and  energy.  The  diffusion  equations  were  uncoupled  by  assuming  a 
Lewis  number  of  unity.  Thus,  the  above  system  of  equations  which 
was  treated  consists  of  a  coupled  set  of  4  non-linear  partial  differen¬ 
tial  equations  of  7th  order  having  split  boundary  conditions. 

For  the  low  Reynolds  number  regime,  all  terms  of  order  1/Reg 
and  larger  were  retained.  The  method  of  separation  of  variables  re¬ 
duced  the  governing  equations  to  a  set  of  5  coupled  non-linear  ordi¬ 
nary  differential  equations  of  order  8  with  an  "unknown”  range  of 
integration.  These  were  then  solved  numerically  utilizing  a  high 
speed  electronic  computer.  Note  that  by  utilizing  a  mass  transfer 


111 


cooling  process,  the  surface  temperature  can  be  maintained  at  suf¬ 
ficiently  low  levels  so  that  heterogeneous  reactions  need  not  be  con¬ 
sidered  explicitly.  The  heat  transfer,  skin  friction  coefficient  and 
the  effectiveness  of  reducing  each  by  mass  transfer  were  determined 
over  the  hypersonic  low  Reynolds  number  range  of  flight  speeds 
from  15,000  ft. /sec,  to  25,000  ft. /sec.  and  an  altitude  range  from 
175,  000  ft.  to  350,  000  ft.  The  new  results  were  compared  with 
earlier  theoretical  studies  in  the  continuum  regime  and  in  the  low  Rey¬ 
nolds  number  regime  (without  mass  transfer),  and  the  differences 
were  interpreted. 


IV 


LIST  OF  SYMBOLS 


c. 

1 


f 

w 

h 


h. 

1 


skin  friction  coefficient 
mass  fraction  of  species  i 

specific  heat  at  constant  pressure  of  the  pure  species 
binary  diffusion  coefficient 
mass  transfer  parameter  (see  eq.  57) 
static  enthalpy 

static  enthalpy  of  the  i  th  species,  including  enthalpy 
of  formation 


Ah£  enthalpy  of  formation  of  the  i  th  species 

i 

k  coefficient  of  thermal  conductivity 

Kq  nose  curvature 

13 

m  mass  transfer  rate  (  =P  V  ) 

W  W  W 

p  static  pressure 

Q  heat  transfer 

Q  energy  flux  vector 

r  cylindrical  radius 

Ro  nose  radius 

13 

T  temperature 

U  x-component  of  velocity 

V  y-component  of  velocity 

V  mass -weighted  average  velocity  of  the  fluid  mixture 
absolute  flow  velocity  of  species  i 


V 


diffusion  velocity  of  species  i 


V. 

1 

flight  speed 

net  mass  rate  of  production  of  species  i  per  unit  volume 
by  chemical  reaction 

X,  y  coordinate  system 

5  viscous  layer  thickness 

5  shock  detachment  distance 

e  density  ratio  across  shock  (=  /P  S  ) 

p  ordinary  viscosity  coefficient 

P  density  of  fluid  mixture 

partial  density  of  species  i 

Z  viscous  stress  tensor 


Subscripts 


s 

immediately  behind  shock 

w 

body  surface 

1,  2 

functions  of  y  alone 

00 

free  stream 

Dimensionless  Groups 


Q 


Le  = 


w 


Pr  = 


C  p 
P 


P  V  R„ 
Re  =  - 

^  Ps 


Stanton  number 


Lewis  number 


Prandtl  number 


Re yno  Id s  numbe  r 


VI 


INTRODUCTION 


For  slender  ballistic  re-entry  vehicles  and  hypersonic  lifting  re¬ 
entry  vehicles  at  high  altitudes,  it  is  mandatory  that  adequate  informa¬ 
tion  be  available  for  prediction  and  control  of  the  exchange  of  mass, 
momentum  and  energy  between  the  surface  of  the  vehicle  and  its  en¬ 
vironment. 

Typical  ballistic  and  lifting  re-entry  trajectories  are  shown  in 
Figure  1.  Note  that  Figure  1  is  typical  for  a  nominal  nose  radius 
equal  to  one  tenth  of  a  foot;  for  smaller  nose  radii,  all  of  the  low 
density  effects  are  shifted  to  lower  altitudes. 

Descriptive  stagnation  region  velocity  profiles  are  shown  to  the 
right  of  the  respective  flight  regimes.  These  regimes  are  discussed 
below. 

As  is  well  known,  hypersonic  flight  at  low  altitudes  is  associated 
with  high  Reynolds  numbers,  and  a  flow  field  defined  by  a  very  thin 
shock  wave  followed  by  an  influscid  region  separated  from  the  body  by 
a  thin  boundary  layer.  This  is  called  the  ordinary  continuum  regime, 
and  the  whole  flow  field  is  known  as  the  shock  layer.  Conditions 
directly  behind  the  shock  wave  are  determined  using  the  Rankine- 
Hugoniot  relations  and  the  influscid  flow  is  determined  by  using  one  of 
the  available  hypersonic  shock  layer  approximations  (see,  for  example 
Ref.  1);  e.g.  modified  Newtonian,  or  stream  tube,  in  order  to  obtain 
outer  boundary  conditions  for  the  boundary  layer. 

For  suborbital  flight  speeds,  the  radiative  heat  transfer  from  the 
shock  layer  is  very  small  compared  with  the  aerodynamic  heat  trans¬ 
fer  and  hence  the  dominant  energy,  momentum  and  mass  transfer  proc¬ 
esses  can  all  be  treated  by  studying  the  phenomena  within  the  boundary 
layer  adjacent  to  the  surface.  This  procedure  has  yielded  many  useful 
-solutions  to  the  thermal  protection  problem,  (e.  g.  Refs.  2-5). 
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At  a  somewhat  higher  altitude  (lower  Reynolds  number)  the  fluscid 
effects  extend  throughout  a  major  portion  of  the  shock  layer  and  the 
flow  may  no  longer  easily  be  separated  into  three  distinct  regions. 
Convention  has  now  termed  this  the  hypersonic  low  Reynolds  number 
flow  regime  which  may  be  divided  into  two  subregimes: 

1)  Viscous  layer  -  It  is  still  assumed  that  the  classical 
conservation  laws  of  macroscopic  flow  are  applicable  and  the  shock 
remains  sufficiently  thin  to  be  treated  as  a  discontinuity  followed  by 
at  least  a  small  influscid  region. 

2)  Merged  viscous  flow  -  At  sufficiently  low  Reynolds  num¬ 
bers  the  now  thickened  shock  merges  with  the  fully  fluscid  flow  behind 
it,  and  the  entire  flow  field  from  the  free  stream  to  the  body  surface 
should  be  treated  as  a  single  problem. 

At  extreme  altitudes,  (the  free  molecule  flow  regime)  it  is  per¬ 
missible  to  neglect  collisions  between  incident  and  reflected  molecules 
in  studying  interactions  between  the  vehicle  and  its  environment.  How¬ 
ever,  it  is  vitally  important  that  the  state  of  the  undisturbed  gas 
(Ref,  6),  and  the  surface-particle  interactions  (Refs,  7-9),  be  fully 
taken  into  account.  With  decreasing  altitude,  particles  reflected  from 
the  surface  begin  interacting  with  the  incident  particles.  It  becomes 
necessary  to  utilize  kinetic  theory,  e,  g,  ,  the  Boltzmann  equation 
(Refs,  10- IZ)  in  order  to  obtain  solutions  in  the  near  free  molecule 
flow  regime. 

Between  the  low  Reynolds  number  and  near  free  molecule  flow 
regimes,  a  transition  regime  exists  which  requires  a  more  rigorous 
kinetic  theory  approach,  possibly  coupled  with  a  macroscopic  analysis 
in  which  higher  order  terms  are  introduced.  It  should  be  noted  that 
the  uncertainty  in  the  extent  of  this  region  is  indicated  by  the  lightly 
shaded  area  plus  the  broken  line  in  the  velocity  profile. 
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It  is  seen  from  Figure  1  that  knowledge  about  the  hypersonic  low 
Reynolds  number  flow  regime  is  important  both  for  slender  ballistic 
entry  vehicles  and  efficient  lifting  re-entry  vehicles.  The  earliest 
studies,  restricted  to  order  of  magnitude  analyses,  are  summarized  by 
Hayes  and  Probstein  (Ref.  1)  where  the  entire  hypersonic  flight  regime 
is  divided  into  seven  subregimes  from  continuiom  to  free  molecule  flow, 
including  boundary  layer  regime,  vorticity  interaction  regime,  viscous 
layer  regime,  incipient  merged  layer  regime,  fully  merged  layer  and 
transitional  layer  regimes,  first  order  collision  theory  regime  and 
free  molecule  flow  regime.  This  classification  scheme  has  been  the 
basis  of  most  subsequent  studies. 

In  Hoshizaki's  study  (Ref.  13)  of  the  effect  of  bluntnes s -induced 
vorticity  on  the  boundary  layer,  it  was  concluded  that  the  inviscid 
momentum  equations  evaluated  behind  the  shock  are  equivalent  to  using 
the  vorticity  as  a  boundary  condition,  and  thus,  the  complete  momen¬ 
tum  equations  incorporate  the  vorticity  boundary  conditions.  Increases 
in  heat  transfer  and  skin  friction  were  obtained  when  the  vorticity 
terms  were  retained  in  the  Navier -Stokes  equations,  as  the  Reynolds 
number  was  decreased.  However,  these  terms  are  retained  in  an 
appropriate  low  Reynolds  number  order  of  magnitude  analysis.  In 
Refs.  14-17,  similar  results  were  obtained  by  appealing  to  the  low 
Reynolds  number  nature  of  the  flow  without  special  recourse  to  vor¬ 
ticity  except  as  a  modification  of  the  pressure  term. 

In  Refs,  13-16,  the  density  was  assumed  constant  through  the 
shock  layer.  Since  it  can  easily  be  shown  that  the  variation  of  density 
through  this  shock  layer  is  of  the  same  order  of  magnitude  as  that 
across  the  shock  wave,  it  would  appear  that  the  assumption  of  constant 
density  is  not  realistic.  In  a  more  recent  study  (Ref.  17),  the  density 
variation  was  taken  into  account,  and  it  was  shown  that  as  the  Reynolds 


3 


number  decreases,  the  ratio  of  heat  transfer  calculated  with  density- 
variation  to  the  heat  transfer  calculated  with  constant  density  increases. 
Clearly,  then,  the  variation  of  density  across  the  viscous  layer  should 
be  considered  in  a  realistic  analysis. 

In  all  of  the  above  studies  however,  solutions  were  obtained  for 
the  flo-w  of  a  viscous  non-dissociating  perfect  gas.  Hence,  these 
studies  ignored  the  important  physicochemical  phenomenon  of  dissocia¬ 
tion,  which  decreases  the  molecular  weight  of  the  gas  and  alters  the 
relationship  between  density,  temperature  and  pressure.  Also  not 
considered  was  the  alteration  in  the  flow  produced  by  mass  transfer 
from  the  surface,  which  will  be  shown  to  have  a  considerable  effect 
upon  the  structure  of  the  viscous  layer  and  other  dependent  variables. 

In  Ref.  14,  Probstein  and  Kemp  extended  their  analysis  to  the 
merged  viscous  layer  problem  by  matching  a  simplified  perfect  gas 
shock  structure  solution  with  their  constant  density  perfect  gas  viscous 
layer  solution.  Recently,  Levinsky  and  Yoshihara,  Ref,  18,  also 
studied  the  merged  viscous  layer  but  included  variable  density. 
Qualitatively,  the  results  appear  credible.  However,  in  addition  to 
the  remarks  made  previously  about  the  non -dissociating  perfect  gas 
model  in  the  viscous  layer  solution,  its  use  in  the  study  of  hypersonic 
shock  wave  structure  in  air  is  even  more  unrealistic,  since  not  only 
is  the  gas  within  a  hypersonic  shock  wave  not  in  chemical  equilibrium, 
it  is  not  even  in  thermodynamic  equilibrium,  (Ref.  19)  i.  e.  ,  the  re¬ 
laxation  times  for  the  translational  degrees  of  freedom  are  much 
shorter  than  the  other  degrees  of  freedom,  (e.  g.  rotational,  vibrational). 
In  particular,  the  available  kinetic  energy  of  the  molecules  is  redis¬ 
tributed  among  many  different  internal  modes  of  energy  storage  which 
absorb  energy  at  different  rates,  but  nevertheless  in  parallel  processes  • 
and  not  in  series.  It  is  stressed  that  the  dissociation  and  ionization 
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processes  which  absorb  energy  and  act  to  reduce  the  gas  temperature, 
will  begin  before  thermal  equilibration  is  complete. 

The  present  analysis  has  been  carried  out  in  the  regime  defined 

in  Figure  1  as  the  hypersonic  low  Reynolds  number  viscous  layer 

regime.  In  this  regime,  the  viscous  layer  thickness  6  is  significant 

compared  to  the  shock  detachment  distance  6  ,  but  6/6  <1,  and 

s  s 

the  shock  wave  thickness  is  sufficiently  small  so  that  it  can  be  as - 
sirmed  that  the  flow  is  in  thermochemical  equilibrium  behind  the  shock, 
before  entering  the  viscous  layer. 

The  governing  equations  utilized  in  the  analysis  are  the  equations 
of  change  for  the  flow  of  a  compressible  chemically  reacting  gas. 
Included  are  the  conservation  of  mass,  momentum  and  energy.  In 
order  to  simplify  the  system  of  equations,  the  diffusion  equations  were 
uncoupled  by  asstrming  a  Lewis  number  of  unity.  The  boundary  con¬ 
ditions  are  split;  some  are  specified  on  the  surface  and  the  others 
behind  the  shock  wave. 

To  evaluate  the  remaining  transport  coefficients,  the  Prandtl  number 
was  assumed  to  be  constant  at  a  value  of  0.  7  1  and  Sutherland's  law  was 
utilized  for  calculating  the  viscosity. 

Even  after  expanding  the  simplified  version  of  the  governing 
equations  into  the  usual  body- oriented  coordinates,  it  is  necessary  to 
make  further  simplifying  assumptions  in  order  to  reduce  the  system 
of  equations  to  a  more  tractable  form.  It  is  further  assumed  that 
the  shock  layer  thickness  §  is  small  compared  with  a  characteristic 
dimension  of  the  vehicle;  in  this  case  the  nose  radius,  R„.  In  fact, 
it  has  been  shown,  e.  g.  Refs.  1,  17,  22, 

^=0(0.1)  (1) 

B 
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in  the  stagnation  region  of  a  hypersonic  vehicle.  Based  on  this  as- 
sumption,  the  order  of  magnitude  of  each  term  in  the  expanded 
governing  equations  may  be  calculated.  Retaining  all  terms  from 
order  one  to  those  of  order  (Rcg)  ^  yields  the  equations  appropriate 
for  a  low  Reynolds  number  analysis.  Following  a  procedure  similar 
to  that  utilized  by  Probstein  and  Kemp  (Ref,  14),  Ho  and  Probstein 
(Ref.  17),  and  Levinsky  and  Yoshihara  (Ref.  18),  the  four  non-linear 
partial  differential  equations  were  reduced  to  an  eighth  order  set  of 
five  non-linear  ordinary  differential  equations  with  split  boundary 
conditions.  The  primary  additional  assumptions  introduced  were: 

1)  the  shock  is  concentric  with  the  body,  2)  the  flow  is  locally  similar 
and  3)  the  variables  are  separable.  The  following  system  of  equations 
represents  the  viscous  flow  in  the  forward  stagnation  region  of  the 
hypersonic  low  Reynolds  number  viscous  layer. 

BASIC  RELATIONS 

The  validity  of  the  conservation  equations  in  the  continuum  regime 
has  been  amply  demonstrated  over  the  past  fifty  years.  In  extending 
the  analysis  to  lower  Reynolds  numbers  it  is  assumed  that  the  density 
is  sufficiently  large  that  the  same  macroscopic  governing  equations 
are  still  applicable.  In  the  absence  of  external  force  fields  the  steady 
state  form  of  these  equations  is: 

Conservation  of  Species  i 

V  '  (p.  V.)  =  W.  (2) 

Conservation  of  Momentum 

p  (v  •  V)  V  =  -  Vp  +  V  .  ^  (3) 
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Conservation  of  Energy 


PV  •  V 

Equation  of  State 


-V-  Q  +  V-  (  r  .  v) 


P  =  P  (h,  p) 


(4) 


(5) 


Summation  over  all  species  yields  the  global  continuity  equation 

V.  (pv)  =  0  (6) 

Implicitly  assumed  In  the  above  definitions  is  that  the  gaseous  system 
is  composed  of  a  mixture  of  chemically  reacting  perfect  gases.  Upon 
neglecting  the  Dufour  effect  and  the  radiative  transport  of  energy,  the 
heat  flux  vector  Q  is  (Ref.  23), 


Q  =  -  kVT  +  ^PiVihi 
i 


(7) 


i.e.  the  sum  of  the  Fourier  conduction  and  the  energy  transferred  by 
diffusion.  From  the  definition  of  h,  it  is  seen  that 


VT  =  =2^ 
C 


Vh-  )'h.  VC 
V  1  1 


where  C  is  the  frozen  specific  heat: 
P 


(8) 


Cp  =  (,) 

As  a  first  approximation,  it  is  assumed  that  the  gaseous  components 
with  similar  thermal  and  chemical  behavior  may  be  treated  as  a  single 
averaged  component  so  that  the  sum  total  can  be  treated  as  a  binary 
mixture  of  "air  atoms"  and  "air  molecules".  Then,  neglecting  pres¬ 
sure  and  thermal  diffusion  effects,  the  mass  flux  with  respect  to  the 
mass-averaged  velocity  reduces  to  the  classical  Fick's  diffusion 


Q  =  -  -tL 

Pr 


Substitution  of  eqs.  (8)  and  (10)  into  (7)  yields 

Vh+  (Le-1) 

Assuming  a  Lewis  number  of  unity  thus  yields 

Q  =  - Vh 

Pr 


(11) 


(12) 


and  the  diffusion  equations  (2)  are  uncoupled,  even  for  the  case  of 
non- equilibrium  flow,  without  resorting  to  any  specific  constraint  on 
the  chemical  source  term,  i,  e,  the  term  Sw  h  is  accounted  for 

i  i  i 

implicitly. 

With  these  simplifying  assumptions  the  governing  equations  appear 
expanded  below  for  the  coordinate  system  shown  in  figure  2. 

Continuity: 

8 


8 

8X 


(pur^)  +  —  (pvr^)  +  K  r^pv=  0 

B 


(13) 


x-Component  of  Momentum: 

8u 


u  +  V 

8x 
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8u  1  8 
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8U 
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y-Component  of  Momentum: 
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8V 

9v  2 

8p  4  8 

- 
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8x 

V  -r—  -  K  U 

9y  B 
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Energy: 
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8x 


h  + 


u 


+  V 


9y 
) 

B  Pr  8y^  '  ^ 


h  + 
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_  8  I  8h_ 

~8y  1 Pr  8y 


(16) 
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It  should  be  noted  that  the  governing  equations  presented  here  are 

more  complete  than  those  given,  for  example,  in  Ref.  17,  in  that  all 

terms  ranging  from  order  unity  to  order  (Re  )  appear.  Although 

s 

Ho  and  Probstein  state  that  they  have  considered  all  terms  ranging 

from  order  unity  to  order  (e  Re  )  actually  one  term  of  order 

-1  ^  . 

(CRe  )  seems  to  be  missing  in  the  energy  equation.  However,  this 

introduces  no  error  in  their  computations  since  the  missing  term  van¬ 
ishes  identically  at  the  stagnation  point. 

State : 

P  =  P  (h,  p)  (17) 


Viscosity  Law: 


11=  11  {h,  p) 


(18) 


At  this  point  then,  there  are  six  unknowns  u,  v,  h,  p,  p  and  p 
which  may  be  determined  by  solving  the  six  equations  (13)  through  (18) 
simultaneously. 


BOUNDARY  CONDITIONS 

The  location  of  the  edge  of  the  viscous  layer  within  the  shock 
layer,  and  the  magnitude  of  the  physical  variables  at  the  outer  edge 
of  the  viscous  layer  are  not  known  "a  priori".  Therefore,  for  con¬ 
venience,  it  was  decided  to  integrate  the  governing  equations  from  the 
wall  to  the  shock  wave,  where  the  physical  variables  can  be  calculated. 

It  was  assumed  that  the  shock  in  the  stagnation  region  is  a  sur¬ 
face  of  discontinuity  concentric  with  the  body.  It  was  further  assumed 
that  thermochemical  equilibrium  is  achieved  behind  the  shock.  Use  of 
the  Rankine -Hugoniot  relations  yield  the  outer  boundary  conditions 
behind  the  shock  wave: 

u  (x,  6  )  =  K  xV  (19) 

s'  B  CO 


9 


(20) 


(21) 

(22) 


wlie  re 


(23) 


It  was  also  assumed  that  there  is  no  slip  or  temperature  jump  at 
the  wall  and  hence  the  boundary  conditions  at  the  wall  are: 


u(x,0)  =11^  =  0 


(24) 


V  (x,  0)  =  V 

w 


(25) 


h  (x ,  0)  =  h 
w 


(26) 


Note  that  the  number  of  boundary  conditions,  namely  seven,  equals  the 
order  of  the  mathematical  system. 


Although  the  magnitude  of  the  physical  variables  can  be  calculated 
behind  the  shock,  the  distance  from  the  wall  to  the  shock  is  not  known. 
However,  one  may  introduce  a  constraint  on  the  family  of  solutions 
which  is  based  on  the  conservation  of  mass  and  takes  the  form: 


V  r 

00  00 


+  P 


V  r  = 
WWW 


6 

r  ®  i 
Jq  pu(2r)''  dy 


(27) 


This  equation  serves  to  delimit  8  . 
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SEPARATION  OF  VARIABLES 


The  governing  partial  differential  equations  may  be  solved  by 
means  of  the  technique  Icnown  as  separation  of  the  variables,  which 
has  previously  proven  useful  in  treating  the  low  Reynolds  number 
viscous  layer,  (Ref.  14).  In  the  manner  utilized  by  Probstein  and 
Kemp,  the  following  separable  forms  of  the  dependent  variables  were 
assumed  which  were  suggested  by  the  boundary  conditions  at  the 
shock  wave. 


u  =  Kg  X  (y) 

(28) 

V=  (l-^Kg^x^)v^(y) 

(29) 

h  =  (1-Kg^  x^)h^(y) 

(30) 

p  =  (1-  ^Kg^x^)  p^(y) 

(31) 

:  (1-Kg^x^)  p^(y)  -  Kg^x^P2(y) 

(32) 

P  =  (1-  Kg^x^)  p^(y) 

(33) 

r  =  X  (1+Kgy) 

(34) 

where  the  terms 
For  convenience, 


with  subscripts  1  or  Z  are  functions  of  y  alone, 
the  following  non-dimensional  forms  are  introduced: 


U  = 


Voo 


- 

v  = 


eVoo 


,  h  = 


__Pl  _  _  ^2 


(35) 


'00"  00 
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where  the  subscripts  co  and  s  refer  respectively  to  the  free  stream 
conditions  and  to  the  state  immediately  behind  the  shock.  Note  that 
the  Reynolds  number  is  defined  by: 


Re 

s 


P  V 
00  00 


(36) 


With  the  above  simplifications  the  low  Reynolds  number  viscous 
layer  equations  reduce  to: 


(  p  v)'  +  (1  +  j)  —  (u+  e  v)  =  0 
£ 


(37) 


U  (u  +  ev)  +  evu' 
pv(ev)'  = 


=  2  (p^  +P2)  +  (p  u')’ +  (2 +i)  /ill'  -  (/iu)'  (38) 


P  (ev)' 


'  +  (lu) 


/i  U 


PJ:l(a4-,v)  +  P2'=  0 


Pr  p  V  h  '  -  Id  h  '  '  +  (1  +  j)  p  h ' 


(39) 


(40) 


(41) 


P  =  p  (h,  p) 


(42) 


p  =  p  (h,  p) 


(43) 


where  the  primes  denote  derivatives  with  respect  to  y.  It  is  seen 
that  separation  of  the  variables  has  increased  the  number  of  unknowns 
to  seven,  including  u,  v,  h,  p^,  p^,  p  ,  jj  •  Similarly,  the  number 
of  governing  equations  is  now  also  seven,  while  the  overall  order  of 
the  mathematical  system  is  eight  requiring  eight  boundary  conditions. 
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The  boundary  conditions  become: 


At  y  =  0 ; 


u  =0,  V  =  v,h  =  h 
w  ’  w  w  w  w 


(44) 


At  y  =  5^ 


“s  '  1’  T  '  Pi  Pz  =  “  (45) 


The  constraint  equation  becomes: 


__ 

(  ®  p  u 

2  (1  +  y) 

e 

_  _ 

dy 


(46) 


Here  it  is  noted  that  for  the  special  case  of  zero  mass  transfer,  i,  e. 

p  V  =0,  examination  of  eq.  (37)  indicates  that  the  constraint 
w 


V  ■  =  0  (47) 

w 

becomes  identical  with  eq.  (46).  This  latter  result  follows  from  the 

fact  that  since  u  and  v  are  both  zero  (for  zero  mass  transfer), 
w  w 

unless^  V  also  vanishes,  one  obtains  the  physically  untenable  result 

W  r  1  j 

of  an  infinite  density  gradient  at  the  wall,  (Ref.  17). 

In  general,  the  transport  properties  depend  on  both  temperature 
and  composition.  However,  since  the  diffusion  equations  have  been 
uncoupled,  the  composition  of  the  gas  is  not  calculated  explicitly,  and 
so  it  is  necessary  to  resort  to  further  approximations.  Constant 
Prandtl  and  Lewis  numbers  have  already  been  introduced  and  hence 
the  thermal  conductivity  and  the  diffusion  coefficients  are  not  required 
explicitly.  Thus,  it  remains  to  evaluate  only  the  viscosity  coefficient. 
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It  is  assumed  that  the  gas  is  in  chemical  equilibrium,  or  more 
explicitly  that  the  state  of  the  gas  can  be  determined  from  a  Mollier 
diagram  (e.  g.  Ref,  24).  That  is,  knowledge  of  two  state  functions 
(e.  g.  pressure  and  enthalpy)  uniquely  determines  the  others  (e,  g. 
temperature,  composition,  density,  etc.  ).  It  is  assumed  here  that 
the  viscosity  coefficient  can  be  approximated  by  Sutherland's  formula 
for  air  (Ref,  25) 


11  =  11  (T)  =  1.16  X  10 


5 

717 

T  ] 

'225+T 

1 

492 
\  / 

3/2 


(48) 


with  T  in  °R,  fd  is  given  in  lb.  /ft.  sec.  However,  since  (from  the 
Mollier  diagram) 


T  =  T  (h,  p) 


(49) 


then  implicitly 


fl  =  fl  (h,  p)  (50) 

and  the  equation  of  state  may  also  be  obtained  from  the  Mollier 
diagram  i.  e. 

P  =  P  (h,  p)  (51) 


Conditions  behind  the  shock  were  determined  from  shock  tables 
which  incorporate  the  1959  ARDC  atmosphere,  (Ref.  26), 

Solutions  were  obtained  on  a  Pace  Electronics  Associates  high¬ 
speed  electronic  analog  computer  in  the  hypersonic  low  Reynolds  num¬ 
ber  flight  regime.  Calculations  were  carried  out  for  a  range  of  flight 
speeds  from  15,000  ft. /sec.  to  25,  000  ft, /sec,  ,  an  altitude  range  from 
175, 000  ft.  to  350,  000  ft.  for  Reynolds  numbers  which  were  varied  between 

25  and  10  ,  and  amass  transfer  parameter,  -0, 4<f  =5  0.  These  results 

w 

were  compared  with  earlier  theoretical  studies  in  the  continuum  regime  and 
in  the  low  Reynolds  number  viscous  layer  regime  (for  zero  mass  transfer). 
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DISCUSSION  OF  RESULTS 


The  non-dimensional  system  of  equations,  (37)  to  (43),  was 
solved  on  a  Pace  Electronics  Associates  analog  computer.  The 
equations  were  integrated  from  the  wall  out  to  the  shock  wave.  Since 
the  boundary  conditions  are  split  it  is  necessary  to  employ  an  iter¬ 
ative  procedure.  Some  typical  solutions  for  the  viscous  shock  layer 
profiles  of  u,  v,  h  and  p  are  shown  in  Figures  3  to  6.  As  in¬ 
dicated  by  previous  studies  the  viscous  effects  extend  further  into 
the  shock  layer  with  decreasing  Reynolds  number. 

After  numerous  calculations  for  a  wide  range  of  altitudes,  flight 
speeds,  Reynolds  numbers  and  wall  temperatures,  it  was  found  that 
all  of  the  data  could  be  correlated  with  several  standard  non-dimen¬ 
sional  ratios  as  functions  of  the  Reynolds  number  alone.  No  dis¬ 
cernible  trends  due  to  the  other  effects  were  noted. 


Introducing  the  Stanton  number 


Q 


w 


p  V  (H  -h  ) 

^00  00  s  w 


(52) 


all  of  the  heat  transfer  data  for  zero  mass  transfer  were  correlated 
as  a  function  of  Reynolds  number,  and  are  shown  in  figure  7,  In¬ 
cluded  in  this  figure  are  also  some  recent  results  of  other  investi¬ 
gators,  Refs,  14,  17,  18,  27  and  28  for  comparison.  Here,  it  is 
seen  that  the  present  analysis  indicates  a  larger  low  Reynolds  number 
effect  than  predicted  by  the  earlier  investigators.  In  order  to  deter¬ 
mine  if  these  effects  were  due  to  the  low  Reynolds  number  nature  of 
the  flow,  i.e,  inclusion  of  all  terms  up  to  order  (Re^)  \  or  to  the 
gas  model  chosen,  new  solutions  were  obtained  for  both  Ho  and 
Probstein's  (Ref,  17),  and  Ferri,  Zakkay  and  Ting's  (Ref.  27) 
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analyses,  utilizing  the  present  gas  model.  These  results  which  are 
based  on  the  use  of  their  differential  equations,  but  our  gas  model, 
are  shown  in  Figure  7  as  curves  4  and  5.  From  these  calculations, 
it  is  seen  that  the  primary  difference  appears  to  be  the  gas  model, 
since  their  original  results  shown  as  curves  2  and  3  of  Figure  7  are 
shifted  upwards  and  agree  very  well  with  the  authors'  results.  It  is 

3 

also  noted  that  for  Reynolds  numbers  less  than  10  ,  the  modified 
results,  (curves  4  and  5)  do  not  increase  with  decreasing  Reynolds 
number  as  rapidly  as  the  present  computations.  This  result  indicates 
the  need  to  retain  more  higher-order  terms  with  decreasing  Reynolds 
number. 

It  may  be  noted  that  in  checking  the  present  system  of  equations, 

including  the  present  gas  model,  at  high  Reynolds  number  (Re  = 

5  ® 

2,  5  X  10  ),  the  heat  transfer  was  found  to  be  approximately  one  per¬ 
cent  different  from  more  exact  digital  calculations  based  on  using  the 
high  Reynolds  number  boundary  layer  equations  and  variable  fluid 
properties,  (Ref,  29), 

The  calculations  in  this  study  were  carried  out  to  Reynolds  num¬ 
bers  below  the  range  in  which  it  is  expected  that  the  viscous  flow 
layer  will  merge  with  the  shock  wave,  in  order  to  determine  trends. 
From  Refs.  17  and  27,  it  appears  that  the  lower  limit  of  the  viscous 

layer  regime  is  approximately  Re  =  100.  Ho  and  Probstein  (Ref, 

s 

17)  extended  their  viscous  layer  results  into  the  merged  viscous  layer 
regime  by  modifying  a  previous  study  (Ref,  14),  (constant  density 
solution).  In  a  more  recent  study  by  Levinsky  and  Yoshihara  (Ref, 

18) ,  the  same  system  of  equations  and  gas  model  employed  by 
Probstein  was  solved  from  the  free  stream  through  the  merged  shock 
and  viscous  layer  to  the  wall.  The  results  are  naturally  similar, 
except  that  paradoxically  they  are  shifted  by  an  order  of  magnitude 
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to  higher  Reynolds  numbers  (See  Figure  7),  with  some  interesting 

consequences.  At  the  higher  Reynolds  numbers  the  curve  almost 

coincides  with  the  authors'  present  results.  Below  a  Reynolds  num- 
3 

ber  of  about  10  the  trend  is  downward  indicating  merged  flow.  An 
extension  of  the  curve  would  appear  to  merge  with  results  from  a 
very  recent  experimental  study  by  Dewey,  Ref.  Z8.  Dewey's  tests 
were  conducted  at  a  nominal  Mach  number  of  5.8  and  stagnation 
temperature  of  226°F  explicitly  ignoring  dissociation. 


It  would  appear  from  the  divergence  of  the  various  theories  and 

2 

experimental  data,  that  the  low  Reynolds  number  regime  Re  <  10 
is  not  yet  entirely  understood. 


With  regard  to  skin  friction,  one  observes  that  the  surface  shear 
stress  is  given  by: 


T  = 
w 


I  8U 


(53) 


and  the  skin  friction  coefficient  by: 


T  = 
w 


(54) 


The  results  of  the  computations  of  all  skin  friction  data  for  zero 
mass  transfer  are  shown  in  Figure  8.  It  is  seen  in  this  figure,  as  in 
the  preceding  one,  that  the  present  study  shows  a  larger  low  Reynolds 
number  effect  than  earlier  studies. 

Here  it  is  seen  that  the  dissociated  gas  model  employed  by  the 
authors  acts  to  increase  the  skin  friction  coefficient  predicted  by 
Refs.  17  and  27,  (see  curves  3  and  4  of  Figure  8).  However,  at  the 
lowest  applicable  Reynolds  numbers  the  modified  calculations  remain 
twenty  percent  below  the  authors'  results.  This  is  actually  not 
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neglected  in 


surprising  since  most  of  the  terms  of  order  (Re^) 
these  earlier  studies  appear  in  the  equation  for  the  x-component  of 
momentum,  which  is  most  influential  in  the  determination  of  the 
shear  stress. 

In  the  presence  of  mass  transfer,  the  energy  transfer  to  the  wall 
is  given  by  (Ref.  29): 

Q  =  (Q  )  + 

w  w . 

m  =  0 
w 

where  the  Q  's  and  (AQ/Aih)  have  opposite  signs  and  m  =  n  v  . 

w  w  w  w 

In  the  high  Reynolds  number  regime  at  large  flight  speeds  Scala  has 

shown  that(AQ/Am)  for  air  injection  could  be  correlated,  independent 

w 

of  pressure  and  geometry,  by 

when  this  "effectiveness  quotient"  is  defined  at  f  =  -0.4  (large  value 

w 

of  the  mass  transfer  parameter).  Figure  9  indicates  that  the  high 
Reynolds  number  approximation  is  quite  good  down  to  Reynolds  num- 
bers  as  low  as  10'^,  for  the  prediction  of  (AQ/Arh)^, 

For  Reynolds  numbers  less  than  10-^,  the  efficiency  of  the  mass 

transfer  process  increases  with  decreasing  Reynolds  numbers.  Note 

that  the  dimensionless  mass  transfer  rate  f  is  defined  as: 

w 

-f 

W 


AQ 

Am  (55) 

/w 
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The  reduction  in  skin  friction  with  mass  transfer  is  similarly- 
given  by  (Ref.  5), 


T 


=0 


w 


A7 

Am 


(58) 


In  Figure  10  the  variation  of  (Ay/Aih)^  with  Reynolds  number 
is  shown.  It  was  found  that  the  new  results  are  in  agreement  with 
Ref.  5  above  Re  >  4x  lO'^.  Here  again,  as  in  Figure  9>  the  effec¬ 
tiveness  of  the  mass  transfer  process  increases  with  decreasing 
Reynolds  number. 


In  the  high  Reynolds  number  regime,  the  shock  detachment  dis¬ 
tance  is  a  function  of  body  shape,  density  and  velocity,  independent 

of  surface  temperature  and  (moderate)  mass  transfer.  However,  at 
3 

Re  <  4x  10  ,  the  shock  detachment  distance  is  a  function  of  body 
s 

shape,  flight  speed,  altitude,  Reynolds  number,  wall  temperature 

and  mass  transfer.  This  effect  was  investigated  in  order  to  predict 

the  effects  of  mass  transfer  in  the  low  Reynolds  number  viscous 

layer  regime  upon  shock  wave  detachment  distance.  Figure  11  shows 

the  variation  of  5  /R  with  the  above  parameters  at  a  single  altitude 

s  B 

and  flight  speed.  Intuitively,  one  expects  that  due  to  the  extended 
interaction  region,  mass  injection  will  push  the  shock  further  away 
from  the  body  with  decreasing  Reynolds  numbers.  As  seen  in  Fig¬ 
ure  11,  this  is  indeed  found  to  be  the  case. 


The  asymptotic  value  of  in  the  high  Reynolds  number 

regime,  was  obtained  from  the  theory  developed  by  Li  and  Geiger, 
(Ref.  Z2). 
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CONCLUSIONS 


1.  Calculations  based  upon  a  dissociating  perfect  gas  model 

predict  that  the  heat  transfer  and  skin  friction  will  increase  faster 

with  decreasing  Reynolds  number,  (for  Reynolds  numbers  less  than 
4 

10  )  than  predicted  by  earlier  theoretical  studies  using  a  non-dis¬ 
sociating  perfect  gas  model. 

2.  Retention  of  all  the  higher  order  terms,  up  to  terms  of 
order  (Re^)  ^  indicates  a  slightly  higher  heat  transfer  rate  at  the  low¬ 
est  applicable  Reynolds  numbers,  but  a  more  significant  increase  in 
skin  friction. 

3.  The  benefits  obtained  with  mass  transfer  in  the  high 
Reynolds  number  regime  (reduction  of  skin  friction  and  heat  transfer) 

have  been  shown  here  to  extend  into  the  low  Reynolds  number  regime. 

*  • 

fact,  the  two  effectiveness  quotients,  (AQ/Am)  and  (A  r/Am) 
which  are  indicative  of  the  efficiency  of  the  mass  transfer  process, 
increase  with  decreasing  Reynolds  number. 
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